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ABSTRACT 

The  hydrodynamic  stability  of  an  underwater  station  was  studied 
using  wind  tunnel  models.   Rigid  model  force  and  moment  coefficients 
were  obtained  for  various  orientations  of  the  station  in  three 
dimensions.   The  unsteady  motion  of  the  station  in  cross  currents 
was  investigated  with  a  dynamically  similar  model.   Color  8  mm 
pictures  and  amplitude- time  traces  of  the  unsteady  motion  were 
obtained. 


fcho 
Moi  i   ornia       9C- 


TABLE   OF   CONTENTS 


PUOLEY  KNOX  LIBRARY 

NAVAL  POSTGRADUATE  SCHUUL 

MONTER^CA  93943-5101 


I.      INTRODUCTION 15 

II.      DISCUSSION 19 

A.  DESCRIPTION   OF   MODELS   AND  EXPERIMENTAL   PROCEDURE    .    .  19 

B.  SCALING  CONSIDERATIONS    22 

III.      RESULTS 30 

A.  RIGID  MODEL  DATA    30 

B.  DYNAMIC  MODEL  RESULTS 31 

IV.   ANALYSIS 33 

A.  STATIC  STABILITY  33 

B.  DYNAMIC  STABILITY 34 

C.  COMMENTS  ON  MUS  CONSTRAINTS 35 

BIBLIOGRAPHY  97 

INITIAL  DISTRIBUTION  LIST   98 

FORM  DD  1473 99 


LIST  OF  TABLES 

NUMBER  PAGE 

I.   Values  of  some  physical  Constants 38 

II.   Model  Summary  .......  39 

III.   Drag  Coefficient  Summary 40 

IV.   Lift  Coefficient  Summary 41 

V.   Pitching  Moment  Coefficient  Summary  42 

VI.   Run  Schedule 43 


LIST  OF  FIGURES 

FIGURE  PAGE 

1.  Mock-up  of  the  Manned  Underwater  Station  .  .  ,  . 50 

2.  Coordinate  axis  system  and  hydrodynamicooef  f  iciehts  ...  51 

3.  Yaw  angle  definition  ............  .  .  52 

4.  Pitch  angle  definition  ..53 

5.  Exterior  dimensions  of  MUS  model 54 

6.  Drawing  of  the  aluminum  spar  used  in  the  rigid  model  .  .  .  55 

7.  Weight  distribution  in  the  dynamic  model  .........  56 

8.  Weight  distribution  in  the  dynamic  model  .57 

9.  Fastening  method  used  for  the  dynamic  model 58 

10.  The  Aerolab  three  component  beam  balance  .........  59 

11  -  44  Photographs  of  the  rigid  model  in  the  Aerolab  32  in. 
x  45  in.  wind  tunnel  showing  support  arrangements  and 
model  orientations 

11.  Model  orientation:   yaw  zero  degrees,  pitch  zero  degrees; 
viewed  from  left  side  of  model;  airflow  left  to  right.  .  .60 

12.  Model  orientation;   yaw  zero  degrees,  pitch  zero  degrees; 
viewed  from  upstream  side   ......  .60 

13.  Model  orientation:   yaw  zero  degrees,  pitch  +30  degrees; 
viewed  from  left  side  of  model;  airflow  left  to  right.  .  .  61 

14.  Model  orientation:   yaw  zero  degrees,  pitch  +30  degrees; 
viewed  from  upstream  side 61 

15.  Model  orientation:   yaw  zero  degrees,  pitch  -30  degrees; 
viewed  from  left  side  of  model;  airflow  left  to  right.  .  .62 

16.  Model  orientation:   yaw  zero  degrees,  pitch  -30  degrees; 
viewed  from  upstream  side   ......  62 

17.  Model  orientation:   yaw  30  degrees,  pitch  zero  degrees; 
viewed  from  left  side  of  model;  airflow  left  to  right.  •  .63 

18.  Model  orientation:   yaw  30  degrees,  pitch  zero  degrees; 
viewed  from  upstream  side 63 


19.  Model  orientation:   yaw  30  degrees,  pitch  +30  degrees; 
viewed  from  left  side  of  model;  airflow  left  to  right  .  .  .  64 

20.  Model  orientation:   yaw  30  degrees,  pitch  +30  degrees; 
viewed  from  upstream  side  ...........  64 

21.  Model  orientation:   yaw  30  degrees,  pitch  -30  degrees; 
viewed  from  left  side  of  model;  airflow  left  to  right  ...  65 

22.  Model  orientation:   yaw  30  degrees,  pitch  -30  degrees; 
viewed  from  upstream  side  .................  65 

23.  Model  orientation:   yaw  60  degrees,  pitch  zero  degrees; 
viewed  from  left  side  of  model;  airflow  left  to  right  ...  66 

24.  Model  orientation:   yaw  60  degrees,  pitch  zero  degrees; 
viewed  from  upstream  side  .  ,  ........66 

25.  Model  orientation:   yaw  60  degrees,  pitch  +30  degrees; 
viewed  from  left  side  of  model;  airflow  left  to  right  ...  67 

26.  Model  orientation:   yaw  60  degrees,  pitch  +30  degrees; 
viewed  from  upstream  side 67 

27.  Model  orientation:   yaw  60  degrees;  pitch  -30  degrees; 
viewed  from  left  side  of  model;  airflow  left  to  right  ...  68 

28.  Model  orientation:   yaw  60  degrees,  pitch  -30  degrees; 
viewed  from  upstream  side  .......  68 

29.  Model  orientation:   yaw  90  degrees,  pitch  zero  degrees; 
viewed  from  left  side  of  model;  airflow  left  to  right  ...  69 

30.  Model  orientation:   yaw  90  degrees,  pitch  zero  degrees; 
viewed  from  upstream  side  ..................  69 

31.  Model  orientation:   yaw  90  degrees,  pitch  +30  degrees; 
viewed  from  left  side  of  model;  airflow  left  to  right  ...  70 

32.  Model  orientation:   yaw  90  degrees,  pitch  +30  degrees; 
viewed  from  upstream  side  ....  f  ............  70 

33.  Model  orientation:   yaw  90  degrees,  pitch  -30  degrees; 
viewed  from  left  side  of  model;  airflow  left  to  right  ...  71 

34.  Model  orientation:   yaw  90  degrees,  pitch  -30  degrees; 
viewed  from  upstream  side  ..................  71 

35.  Model  orientation:   yaw  zero  degrees,  pitch  +60  degrees; 
viewed  from  left  side  of  model;  airflow  left  to  right  ...  72 

36.  Model  orientation:   yaw  zero  degrees,  pitch  +60  degrees; 
viewed  from  upstream  side  . 72 


37.  Model  orientation:   yaw  zero  degrees,  pitch  +90  degrees 
viewed  from  left  side  of  model;  airflow  left  to  right  ...  73 

38.  Model  orientation:   yaw  zero  degrees,  pitch  +90  degrees; 
viewed  from  upstream  side  ........  73 

39.  Model  orientation:   yaw  zero  degrees,  pitch  -60  degrees; 
viewed  from  left  side  of  model;  airflow  left  to  right  ...  74 

40.  Model  orientation:   yaw  zero  degrees,  pitch  -60  degrees; 
viewed  from  upstream  side 74 

41.  Model  orientation:   yaw  zero  degrees,  pitch  -90  degrees; 
viewed  from  left  side  of  model;  airflow  left  to  right  ...  75 

42.  Model  orientation:   yaw  zero  degrees,  pitch  -90  degrees; 
viewed  from  upstream  side  ........  75 

43.  Support  arrangement  used  for  yaw  variation  in  five  degree 
increments;  model  at  zero  degrees  pitch;  viewed  from 

side;  flow  left  to  right ...76 

44.  Same  as  43,  but  view  from  upstream  side 76 

45.  Drag  coefficient  contour  plot  ,..,...,  77 

46.  Drag  coefficient  vs  Reynolds  number , 78 

47.  Pitching  moment  coefficient  vs  Reynolds  number  .......  79 

48.  Lift  coefficient  vs  Reynolds  number 80 

49.  Drag  coefficient  vs  pitch  angle;  yaw  zero  degrees 81 

50.  Pitching  moment  coefficient  vs  pitch  angle;  yaw  zero 
degrees.   Reynolds  number  held  constant  .  .  82 

51.  Lift  coefficient  vs  pitch  angle;  yaw  zero  degrees. 

Reynolds  number  held  constant .  .  .  83 

52.  Conversion  scale.  ................  84 

53.  Yawing  moment  coefficient  vs  yaw  angle;  pitch  zero 

degrees.   Reynolds  number  held  constant  ..........  85 

54.  Drag  coefficient  vs  yaw  angle;  pitch  zero  degrees. 

Reynolds  number  held  constant 86 

55.  Side  force  coefficient  vs  yaw  angle;  pitch  zero 

degrees.   Reynolds  number  held  constant  ...  87 

56.  Photograph  of  the  West  Coast  Research  48  in.  x  60  in.  wind 
tunnel,  showing  support  beams  for  the  wire  attachment 
exterior  to  the  test  section   ..........88 


57.  The  dynamic  model  in  the  test  section,  side  view  .....   88 

58.  Photograph  of  the  dynamic  model  in  the  test  section, 

view  from  downstream  , .  ..........   89 

59.  Close-up  of  the  dynamic  model   ..............   89 

60.  Photograph  of  dynamic  model  instrumentation  hardward  and 
bearings  .  .  ...........   90 

61.  Photograph  of  the  modified  potentiometer,  used  in  dynamic 
model  instrumentation   .........  ,   90 

62.  Damped  frequencies  vs.  dynamic  pressure,  for  model  C  in 

yaw  mode 91 

63.  Damped  frequencies  vs.  dynamic  pressure,  for  model  D  in 

yaw  mode  ..................  92 

64.  Damped  frequencies  vs  dynamic  pressure,  for  model  E  in 

yaw  mode 93 

65.  Damped  frequencies  squared  vs.  dynamic  pressure  for  the 

yaw  mode 94 

66.  Damping  ratio  vs  dynamic  pressure  for  the  yaw  mode  ....   95 

67.  Time  history  of  free  vibration  in  yaw,  strip  recorder 

traces .  .   96 


10 


LIST  OF  SYMBOLS 

k  Projected  frontal  area  (characteristic  area) 

oC  Pitch  angle 

c<  Pitch  angle  rate,  ^  °S 

dt 
Force  Coefficients: 

Cp  Drag  coefficient,  _EL 

^cr  Side  force  coefficient,   £_ 

C  Lift  coefficient,   L._ 

L  %A 

Moment  Coefficients: 

CM  Pitching  moment  coefficient,  i_L 

C^  Yawing  moment  coefficient,  U 

Stability  derivatives:  ° 


CN\  .  ^  *~hA 


H 


be* 


KWi 


CF„.  3CF 


J" 


Vf> 


cr-.  *c? 


11 


D  Drag  force 

ci  Mean  diameter  (characteristic  length) 

<  Side  force 

9r  Gravitational  constant 

1^  Mass  moment  of  inertia 

-$?.  Spring  constant 

A  Scale  ratio,  _moJel 


K 


c^-?ro-\o4t1f>€ 


1—  Lift   force 


Pitching  moment 


~*~n-  Mass 

'V  Yawing  moment 

t\/Re  Reynolds  number,  Y_d 

~L>  Kinematic  viscosity 

— ^  Dimensionless  kinematic  viscosijy, 

^salt  water 


~0  air 
5  Damping  ratioy.  dimensionless 

P  Density 

-i-Pv2" 

Si-  Dynamic  pressure,  2. 

~Y  Tension  of  cable 


£  Time 

x  /  ZV 


*£*  Dimensionless  time,,  L  {,  "3T 


V  Velocity 

V  Dimensionless  velocity, 


V 


model 


Vprototype 
W  Weight 

£-0  Natural  frequency  of  yaw  mode 


12 


£-0^  Dimension  less  frequency  of  yaw  mode  =  iL  ^r\ 

Y  Yaw  angle 

"1^  Yaw  angle  rate,  ^1^ 

X  V  2f  Coordinates,  displacement 

-Tl^  Natural  frequency  of  lateral  translational 

mode 
— CL^  Dimensionless  frequency  of  lateral 

translational  mode 
Subscript : 
yrL  model 

D  prototype 


13 


ACKNOWLEDGMENT 

The  writer  wishes  to  express  his  sincere  appreciation  to  those 
persons  whose  assistance  and  cooperation  contributed  greatly  to  the 
completion  of  this  study,  in  particular,  to  Mr.  Gilbert  DeVlaminck, 
who  fabricated  the  MUS  models,  to  Mr.  Stan  Johnson  and  Mr.  Ted  Dunton 
of  the  Naval  Postgraduate  School  Department  of  Aeronautics,  and  to 
the  technicians  of  the  Department  of  Aeronautics. 

The  writer  is  especially  indebted  to  his  thesis  advisor, 
Dr.  Louis  V.  Schmidt,  for  his  guidance  in  this  research  effort. 


14 


I,   INTRODUCTION 

The  purpose  of  this  study  was  to  investigate  the  hydrouynamic 
stability  of  a  manned  underwater  station  when  the  station  is  fully 
submerged  and  subject  to  cross  currents.   The  concept  of  a  manned 
underwater  station,  or  MUS,  has  been  investigated  by  General  Dynamics 
Corporation,  Electric  Boat  Division,  under  contract  for  the  U.  S. 
Naval  Civil  Engineering  Laboratory,  Port  Hueneme,  California  (Ref.  1). 
The  present  study,  which  is  in  support  of  the  earlier  MUS  work,  was 
sponsored  by  NCEL  under  work  order  WR-8-0128  and  was  monitored  by 
Dr.  Cheng  Lung  Liu  of  NCEL  Ocean  Engineering  Department- 

The  object  of  the  MUS  project  as  stated  in  Ref.  1  is  to  develop 
a  deep  ocean  submergence  habitat  which  will  accommodate  a  five  man 
crew  for  an  operational  period  of  thirty  days.   The  station  is  to 
be  capable  of  submergence  to  depths  of  6,000  feet. 

As  envisioned  by  General  Dynamics,  the  station  would  consist  of 
two  circular  cylinders  of  unequal  diameters  and  heights.   The  smaller 
diameter  cylinder  is  capped  with  hemispherical  ends  and  will  contain 
a  nuclear  reactor.   The  larger  diameter  cylinder  is  capped  with 
spherical  ends  which  are  fitted  to  the  cylinder  by  conical,  tangential 
shells.  (Fig.  1)   The  habitat  cylinder  will  have  five  decks  and  is 
the  working  space  for  the  crew. 

Since  these  cylinders  are  designed  only  for  a  vertical  orientation, 
a  large  ballast  tank  completely  encircles  both  cylinders.   The  design 
of  this  tank  is  based  in  part  on  weight  and  balance  considerations. 
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The  station  is  not  equipped  with  a  propulsion  device  although 
it  does  have  ballast  and  buoyancy  control  equipment.   The  concept 
of  the  station  is  somewhat  unique  and  the  MUS  may  best  be  described 
as  an  elevator.   In  concept,  the  station  is  to  be  towed  on  the  surface 
to  a  pre-selected  site  by  a  support  vessel.  When  the  operations 
site  is  reached,  the  station  will  be  manned  and  an  emplacement 
procedure  will  commence.   The  MUS  will  utilize  an  anchor  and  cable 
arrangement  to  pull  itself  to  the  bottom  in  bootstrap  fashion.   It 
is  equipped  with  a  winch  (located  on  the  main  ballast  tank)  and 
attendant  systems  for  control. 

When  submerged,  the  MUS  will  maintain  a  slight  positive 
buoyancy  at  all  times.   The  anchor  will  first  be  lowered  to  the 
ocean  bottom  while  the  MUS  holds  on  the  surface.   Once  the  anchor 
is  in  place,  the  buoyancy  of  the  MUS  is  reduced  and  winchdown  begins. 
As  the  winch  retrieves  the  cable,  the  MUS  will  haul  itself  to  the 
bottom.   The  recovery  cycle,  i.e.,  the  ascent  to  the  surface,  is 
just  the  reverse  of  this  procedure. 

It  is  primarily  with  the  emplacement  and  recovery  operations  that 
this  study  is  concerned.   When  the  MUS  is  in  the  mid-water  column, 
it  may  well  be  subjected  to  a  cross  flow  caused  by  deep  ocean 
currents.   Evidence  exists  to  expect  such  currents  throughout  the 
water  column,  having  a  magnitude  of  0,1  to  0.5  knots, 

If  such  currents  were  to  exist,  it  would  be  reasonable  to  expect 
the  MUS  to  be  forced  into  motion  until  equilibrium  is  attained. 
Such  motion  could  be  translational  or  rotational  and  perhaps 
oscillatory  in  nature.   Since  the  MUS  is  restrained  only  by  a  single 
cable  (to  the  anchor),  the  only  other  stabilizing  elements  are: 
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1.  The  buoyancy  couple  caused  by  the  displacement  of  the  meta 
center  relative  to  the  center  of  gravity. 

2.  Hydrodynamic  stabilization,  i.e.,  the  geometry  of  the  station, 
including  such  extremities  as  the  legs  and  pads. 

3.  The  inertial  resistance  to  motion  (moment  of  inertia) 
Another  aspect  of  the  MUS  recovery  cycle  which  necessitated  a 

stability  investigation  was  that  of  emergency  ascent.   It  was  proposed 
that,  in  the  event  of  an  emergency  such  as  cable  jamming  or  bottom 
entrapment  of  the  support  pads,  a  free  ascent  could  be  accomplished 
by  cutting  the  cable  (with  explosive  cutters)  and/or  jettisoning  the 
support  pads  and  legs.   Since  these  are  in  fact  stabilizing  devices, 
such  an  ascent  could  only  be  feasible  if  the  resulting  configuration 
of  the  MUS  provided  sufficient  stabilization. 

The  objectives  of  the  experiments  reported  herein  were  twofold, 
name ly : 

a.  To  establish  by  means  of  rigid  model  wind-tunnel  tests  the 
pertinent  dimensionless  force  and  moment  coefficients  for 
the  MUS  at  various  orientations,  and 

b.  To  determine  the  dynamic  behavior  using  a  "partial  mode" 
dynamic  model,  and  in  particular  to  note  the  existence  of 
any  self-induced  and/or  vortex-induced  oscillations. 

In  this  investigation,  the  MUS  was  represented  by  an  0.0209 
scale  model.   Reynolds  numbers  were  matched  between  model  and  prototype 
so  that  aerodynamic  measurements  on  the  scale  model  were  applicable 
to  the  hydrodynamic  traits  of  the  prototype. 

One  should  note  that  the  MUS  configuration  is  still  in  the  early 
design  stages,  hence  the  results  obtained  in  this  study  should  have 
as  a  goal  the  possible  application  to  other  configurations.   However, 
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accurate  forecasting  tends  to  be  uncertain  because  the  bluff-body 
features  of  these  configurations  do  not  interact  in  a  linear  manner. 
The  testing  accomplished  in  this  study  was  relatively  easy,  and  other 
configurations  could  be  studied  using  a  similar  approach. 
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II  DISCUSSION 

A.   DESCRIPTION  OF  MODELS  AND  EXPERIMENTAL  PROCEDURE 

The  investigations  described  herein  upon  the  0.0209  scale  MUS 
models  were  performed  at  the  Naval  Postgraduate  School  (NPS)  low- 
speed  wind  tunnel  facilities.   The  force  and  moment  measurements 
were  conducted  on  the  "rigid"  model  in  the  Aerolab  32  in.  x  45  in. 
wind  tunnel  while  the  scaled  dynamic  model  was  tested  in  the  West 
Coast  Research  48  in.  x  60  in.  wind  tunnel. 

The  "rigid"  model  was  linearly  scaled  and  was  fabricated  from  an 
overlay  of  Honduras  mahoganyy  upon  an  aluminum  spar  or  reinforcing 
plate  (Fig.  6).   The  purpose  of  the  spar  plate  was  to  make  possible 
the  attachment  of  the  support  arms  and  tail  sting  at  prelocated  or 
orientations  of  the  model.   It  was  possible  to  position  the  model 
in  a  multitude  of  ways  upon  the  conventional  three  strut  support 
system  so  that  representative  yaw  and  pitch  attitudes  could  be 
obtained  for  an  angle  range  of  -180  and  -90   respectively.   Figure 
5  presents  a  sketch  of  the  model  including  pertinent  dimensional 
data  relative  to  the  force  model,  while  Figures  11  to  44  show  the 
model  installed  in  the  wind  tunnel  at  various  orientations.   The 
model  was  tested  at  dynamic  pressures  from  3.85  to  57.8  psf  (in  air) 
which  correspond  to  prototype  velocities  in  sea  water  from  0.08  to 
0.30  knots  (refer  to  Figure  52  for  conversion  from  model  to  prototype 
operating  conditions). 

The  forces  and  moments  were  measured  upon  an  Aerolab  three- 
component  beam  balance  (Fig.  10)  which  had  the  capability  of  measuring 
lift,  drag,  and  pitching  moment  relative  to  the  wind  axes. 
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The  balance  limits  were: 

Component  Limits  Sensitivity 

Lift  force  -  150  lbs  -  0.010  lbs 

Drag  force  -  50  lbs  -  0.005  lbs 

Pitch  moment  -  1000  in-lbs  -  0.050  in-lbs 

Care  should  be  exercised  in  interpreting  the  meaning  of  lift  and 
drag  forces  and  pitching  moments  at  some  of  the  extreme  model  positions, 
such  as  when  the  model  has  been  rolled  90  degrees  on  the  balance, 
since  in  this  attitude  pitch  angle  changes  actually  corresponded  to 
model  yaw  angle  changes  (Figs.  43  and  44).   In  this  extreme  case, 
lift  force  and  pitching  moment  correspond  to  side  force  and  yawing 
moment  respectively.   It  is  unfortunate  that  the  balance  system  was 
not  readily  usable  as  a  six-component  type  since  then  it  would  have 
been  possible  to  convert  the  wind  axes  results  to  stability  or  body 
axes  as  is  frequently  done  in  aircraft  analysis. 

Data  corrections  were  made  for  weight  tares,  while  aerodynamic 
tares  to  account  for  the  influence  of  the  support  system  were 
approximately  estimated  by  measuring  the  balance  force  and  moments 
without  a  model  present.   Since  the  model  was  a  relatively  bluff  body, 
as  compared  to  a  streamlined  aircraft  configuration,  these 
approximations  were  considered  as  being  reasonable. 

The  dynamic  model  tests  may  be  described  as  "partial  mode"  tests 
since  only  selected  degrees  of  freedom  were  investigated.   The  modes 
investigated  included  lateral  and  axial  translation  and  yaw  rotation. 
The  tests  were  performed  upon  two  0.0209  scale  models,  one  made  from 
balsa  and  the  other  from  sugar  pine  (Figs,  7,  8,  and  9). 
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The  balsa  model  had  as  its  purpose  an  initial  observation  of  the 
model's  dynamic  behavior  and  an  experimental  confirmation  of  the 
yawing  moment  static  trim  points.   The  full  scale  weight  and  weight 
distribution  table  for  the  MUS  from  Ref.  1  was  used  as  a  basis  for 
design  of  the  pine  model.   By  using  a  component  approach,  the  full 
scale  weights  and  weight  distribution  were  scaled  to  obtain  the  desired 
mass  moment  of  inertia  for  the  pine  model.   Removable  lead  inserts 
were  utilized  to  allow  a  variation  of  mass  and  mass  moment  of  inertia. 
The  balance  of  the  model  about  the  desired  center  of  gravity  location 
was  checked  experimentally,  while  the  mass  moment  of  inertia  was 
determined  by  utilizing  the  principle  of  the  torsional  pendulum. 
An  object  having  a  known  moment  of  inertia,  namely  a  solid  steel 
cylinder  of  given  weight  and  geometry,  was  suspended  on  an  0.037 
inch  diameter  steel  piano  wire  of  known  length.   The  oscillation 
frequency  was  then  measured  which  in  turn  allowed  the  determination 
of  the  wire's  shear  modulus  of  rigidity.   From  the  physical  properties 
of  the  wire  (diameter,  length,  and  shear  modulus),  and  the  measured 
free  vibration  oscillation  frequencies,  it  was  possible  to  determine 
the  mass  moment  of  inertia  of  the  model.   Table  II  lists  the  physical 
properties  of  the  dynamic  models  tested. 

A  narrow  diameter  tube  was  inserted  along  the  vertical  axis 
(through  the  assumed  prototype  center  of  gravity)  in  order  to  shield 
the  dynamic  model  from  the  support  wire.   The  model  was  suspended  on 
an  0.037  inch  diameter  steel  piano  wire  which  was  fixed  to  supports 
external  to  the  tunnel.   (Figs.  56  through  59)   The  wire  was  preloaded 
to  50  lbs.  tension  by  means  of  turnbuckles,  and  the  tension  was 
frequently  checked  by  using  a  standard  rigging  tensiometer.   Initially, 
only  the  taut  wire  support  was  used,  but  later  tests  were  conducted 
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with  the  translational  degrees  of  freedom  removed  by  means  of  local 
tie-down  constraints.   In  addition,  the  last  sequence  of  dynamic 
model  tests  were  provided  with  an  electrical  angle  sensing  potentiometer 
in  order  to  allow  time-history  traces  to  be  made  of  the  angular  motion. 
(Figs.  60  and  61) 

Motion  picture  photography  of  the  dynamic  model  tests  was  accom- 
plished by  using  a  mirror  arrangement.   An  adjustable  mirror  was 
mounted  to  the  tunnel  wall  above  the  model  so  that  overhead  photography 
was  possible.   The  later  tests  involved  simultaneous  photography 
and  electrical  data  collecting.   (Table  VI) 

B.   SCALING  CONSIDERATIONS 

The  scale  of  the  models  (0.0209)  was  established  on  the  basis  of 
the  largest  model  compatible  with  the  wind  tunnel  test  section  size 
and  velocity  range  in  combination  with  a  Reynolds  number  match 
between  model  and  prototype. 

For  the  prototype  design  with  an  ocean  current  of  0..3  knots  in 
salt  water,  an  estimate  of  full  scale  Reynolds  number  would  be 

N/Re=  VpAp  -  (0.1   Krs  x  Uj  Vr")  (  »*Vft3  ^O 
-V?  (oin  »  ,o~4  *^ 

W^g:  =  *f.  4-5"  X  10  (prototype  in  salt  water) 

The  above  Reynolds  number  is  based  upon  a  mean  diameter  of  14.93 
feet.   An  0.0209  scale  model  (  /N  =  0,0209  =  scale  ratio)  would  satisfy 
the  prototype  Reynolds  number  condition  above  at  a  tunnel  velocity  of 

V^=  ^^_K/Rg  *Q.S9   x  /o"4  -Sec- X*'^  *  '<>*)   s  %VSA    Jj* 

cir*         (0.3/3  f+^ 
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■which  corresponds  to  a  tunnel  dynamic  pressure  of  60.3  psf.   The  ratio 
of  model  to  prototype  velocities  is  443  under  the  assumption  of 
equivalent  Reynolds  numbers. 

The  mass  moment  of  inertia  scaling  of  the  dynamic  model  was 
established  by  the  free  vibration  equation  in  the  yaw  mode.   As 
an  approximation,  one  may  consider  the  slope  of  the  C^  vs  y      curve 
to  be  linear  in  the  immediate  neighborhood  of  a  zero  moment  static 
trim  position. 

The  total  yawing  moment  may  be  approximated  by  the  sum  of  the 
linearized  static  term  and  a  term  due  to  the  dimensionless  angle 

rate  (.^V 


cN-    LC»JY  *  h  ^»v~\V 


f  J  t  zv     L       y  A    '  (1) 

Multiplying  ^  by  (ft^d)  to  obtain  physical  units,  then  the  free 
vibration  equation  in  yaw  is: 


dt  r  etc.  7 

(2) 


where:        ^V»-  =  ^B-Ss   moment  of  inertia,  taken  about  the  £  axis. 

;  tl 
where:  £    —        "t 

<&> 

then 

(3a) 


Dividing  through  by  T^       and  introducing  the  dimensionless  time,"£* 
2.V 


a  •'dT1      x&    y-dt-      >  % 


The  enforcement  of  Reynolds  number  equivalence  insures  that  the 
dimensionless  yawing  moment  coefficient  ^~ts/-i         i-s  identical  for 
both  model  and  prototype.   In  order  to  achieve  full  dynamic  similarity, 
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the  nondimensional ,  natural  frequency  of  the  model  must  be  identical 
with  that  of  the  prototype.   This  argument  then  defines  the  scaling 
law  for  the  mass  moment  of  inertia  as  shown  below. 


This  may  be  compared  to  the  familiar  vibration  equation 
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(3b) 
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In  order  for  the  equation  of  motion  to  describe  concurrently  the 
behavior  of  both  model  and  prototype,  each  term  must  equate  in 
model  to  prototype  conversion,  therefore 
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if  the  inertial  scaling  of  above  were  used,  then  the  non-dimensional 
frequency  would  be  alike  between  model  and  prototype,  i.e., 
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Similarly,  a  check  of  the  damping  terms 
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will  yield  the  same  result  by  the  same  argument. 

Unfortunately,  such  scaling  was  not  physically  realizable  as  it 
would  demand  a  prohibitively  light  construction  of  the  model. 

Returning  to  the  non-dimensional  natural  frequency,  which  may 
be  identified  in  eqn.  (3b); 
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if  the  mass  moments  of  inertia  were  set  as 


then  the  dimensionless  frequencies  won't  correspond,  but  will  scale 
as:    /  —  "2- 
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Considering  that  dimensionless  frequency  is  by  definition: 
CO, 
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one  readily  obtains  that 
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After  substitution  of  the  scaling  parameters  (cf.  Table  I),  the  true 
frequency  scaling  becomes: 

(60 ^p  =   O.   OOI34  C^V 

For  example,  if    {^\y^     "  °'  7°   °PS  ,  then 

C^OP  =  0,0013^  (0.70  ^ps^>(60  ^  =  s,oST£2   ^$ 

and  the  prototype  period  would  be  18.1  min/cycle. 

Using  a  quasi -steady  argument,  enforced  by  the  case  where  the 
damping  ratio   F   is  zero,  one  may  use  the  above  scaling  law 
to  convert  model  frequencies  to  prototype  frequencies  for  the  yaw  mode, 

In  a  similar  fashion  to  the  above  analysis  concerning  the  yaw 
mode,  one  may  consider  the  trans lational  mode  free  vibration  equation. 
In  the  lateral  translational  mode,  a  similar  analysis  to  the  yaw 
mode  scaling  may  be  made  starting  from  the  free  vibration  equation, 

3  °  (7) 

The  damping  derivative,   Cp  .    ,  may  be  approximated  by  the  quasi- 

0 
steady  assumption  in  a  manner  similar  to  that  used  by  Parkinson, 

(Ref.  2)  when  he  considered  the  limit-cycle  solution  for  the  galloping 

transmission  line  problem.   Using  the  sign  convention  of  Figure  2, 

one  obtains  that:  I   ^ 

S  =   +  7      V 

« 
since    T^  =   +  3/, 

Substitution  of  this  approximation  into  equation  (7)  yields  that: 

6  r^y   F^  0      mtf  (8) 
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It  is  convenient  to  convert  equation  (8)  to  dimensionless  time,  L      > 
and  dimensionless  lateral  translation,  L\       ,  where   H  —   VA 
Direct  substitution  plus  a  slight  amount  of  algebraic  simplification 


yields 
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Equation  (9)  may  be  compared  with  the  second  order  differential 
equation;  equation  (4),  in  order  to  identify  that  in  the  lateral 
translational  mode,  the  non-dimensional  frequency  is: 


where    Jl   »  /^ 


For   the  model    installation,    the   spring  constant,    k,    is   proportional 
to   the   tension        /        in   taut  wire   support   as   may  be.  seen   in  the   rear 
view  sketch  below: 
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An  approximation  for  the  tethered  prototype  vehicle  may  be  made  on 
the  assumption  that  the  cable  tension  is  equivalent  to  the  buoyancy 
force.   For  this  case,  one  would  obtain 


-ft   =   2  TP 
*p       Xp 
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where  the  tether  cable  length,   JL  p   ,  is  an  equivalent  value  which 
may  be  estimated. 

The  initial  consideration  will  be  to  find  the  true  frequency  value 
of  the  prototype  on  the  assumption  of  equivalent  non-dimensional 
frequencies,  followed  by  a  correction  for  differences  in  cable  tension. 
For  the  case  of   /  —  >.         /  ~  \ 


(^z(Ji^B^ 


(-*«>  =  ih^r  4*"V 


Next,  we  shall  investigate  the  implications  of  cable  tension.   Since 
-it^   ""   ^y.    3  it  follows  that: 

/JfepN*-  =  fine  \*  (  — — N 

^m)      ~~      ^75"    J      ^.  114*13 J  67.8 

.%,     a?  _i_  jfe' 

P  44oo         * 

If  cable  lengbbs  are  approximately  given  by     >C  n   ==  600  ft., 
J?   -  6  ft.  and  model  cable  tension  is  50  lbs.,  then 
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Actually  it  is  envisioned  that  the  prototype  buoyancy  force  will  be 
on  the  order  of  100  lbs.,  hence  the  corrected  frequency  would  be 

Corr  ' 


C-^„    «•   sfe  ^ 


w> 


Therefore  if  the  lateral  translational  frequency  of  the  model  were 
about  1  cps ,  then  the  representative  prototype  period  would  be 
about  53  minutes  per  cycle. 

In  retrospect,  the  dynamic  scaling  problem  would  require 
simultaneous  scaling  of  mass,  mass  moment  of  inertia  about  the 
yaw  axis,  and  cable  tension.   Since  this  procedure  was  not  initially 
followed,  estimates  were  made  to  account  for  the  relative  effects 
of  the  mismatches. 
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III.   RESULTS 

A.   RIGID  MODEL  DATA 

Various  rigid  model  plots  are  included  (Figs.  45  through  55). 
Test  runs  1  through  37  indicated  that  the  hydrodynamic  coefficients 
are  essentially  constant  throughout  the  Reynolds  number  range  under 
consideration.   A  large  number  of  Coefficient-Reynolds  number  plots 
were  obtained  (see  run  schedule,  Table  VI),  and  representative 
plots  are  shown  in  Figures  46  through  48  inclusive.   The  largest 
values  of  the  hydrodynamic  coefficients  that  were  measured,  with  the 
corresponding  full  scale  forces  and  moments,  are  plotted  for 
comparative  purposes  in  Figure  52.   These  forces,  as  compared  to 
the  buoyancy  force  of  100  lbs.,  are  considered  to  be  significant. 

A  drag  coefficient  contour  plot  is  included  in  Figure  45.   The 
drag  coefficient  values  along  the  ordinate  and  absisca  of  this  plot 
were  measured  in  five-degree  increments  (See  also  Figs.  49  and  54), 
while  other  values  were  measured  in  combinations  of  thirty-degree 
increments . 

The  plot  of  yawing  moment  coefficient  versus  yaw  angle  (Fig.  53) 
shows  that  the  MUS  has  two  stable  trim  points  in  yaw,  at  zero  degrees 
and  180  degrees,  while  two  unstable  trim  points  exist  at  90  and  270 
degrees  of  yaw.   It  should  be  emphasized  that  these  are  zero  yawing 
moment  positions,  while  the  slope  indicates  the  degree  of  stability. 
Also  note  that  stable  trim  did  not  occur  at  the  orientation  for 
minimum  drag. 

Direct  measurement  of  the  yawing  moment  vs  yaw  angle  has 
application  to  defining  the  aerodynamic  restoring  moment  term  in  the 

30 


yaw-degree-of-freedom  dynamic  equation.   However,  means  do  not  exist 
for  direct  measurement  of  the  damping  term  due  to  yaw  rate.   The 
damping  term  must  be  deduced  indirectly  by  observing  the  time  history 
of  decaying  yaw  oscillation. 

The  plot  of  side  force  versus  yaw  angle  at  zero  pitch  angle 
is  shown  on  Figure  55.   An  approximate  analytic  expression  for  the 
variation  of  side  force  coefficient  is: 

CF       &      o."36  Hh    f 

However,  the  actual  plot  shows  many  irregularities  in  the  curve 
which  may  be  related  to  the  complex  flow  interactions  between  the 
various  bluff -body  components  that  comprise  the  MUS.   Especially 
noteworthy  is  the  rapid  unstable  trim  points,  "W   =  -90  . 

It  should  also  be  noted  that  the  above  expression  for  side  force 
coefficient  will  result  in  the  damping  term  of  the  lateral  trans- 
lational  degree  of  freedom  having  a  negative  value.   The  consequence 
of  this  feature  was  noted  in  the  dynamic  model  as  a  limit  cycle 
type  of  motion. 
B.   DYNAMIC  MODEL  PxESULTS 

A  color  8  mm  motion  picture  film  was  prepared  as  a  supplement 
to  this  report.   Part  one  of  the  film  shows  the  experiments  with 
the  single  wire  arrangement,  and  careful  observation  will  show  the 
lateral  motion  beating  phenomenon.   In  this  configuration  the  model 
had  degrees  of  freedom  both  in  yaw  and  lateral  translation.   The 
combination  of  an  apparently  lightly  damped  yaw  mode  plus  a 
negatively  damped  lateral  translation  mode  produced  a  complex 
response  with  an  apparent  indication  of  beating  due  to  the  closeness 
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of  the  two  frequencies.   The  lateral  excitation  caused  the  self- 
induced  nature  of  the  oscillation  in  yaw,  but  this  did  not  appear 
when  the  lateral  constraint  was  imposed. 

Part  two  of  the  film  shows  the  instrumented  and  constrained 
model  experiments.   No  lateral  motion  occurred  and  the  damping 
ratio  was  small  but  finite  until  a  dynamic  pressure  of  five  or  six 
psf  was  reached.   Strip  recorder  charts  are  included.   (Fig.  67) 
Dynamic  model  data  are  summarized  below  in  Figs.  62  through  66. 
The  predicted  frequencies  for  the  model,  shown  as  estimates  on 
Figs.  62  to  64,  are  based  upon  the  recognition  of  frequency  from 
equations  (2)  and  (4a)  as: 


?  V      — 


(.%M}  /  c     \ 
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Although  the  predicted  frequencies  are  about  10  to  15  percent  higher 
than  measured  values,  the  trends  with  dynamic  pressure  are  proper 
and  add  confirmation  to  the  concept  of  an  equivalent  hydrodynamic 
torsional  spring  term.   It  is  quite  possible  that  the  actual  observed 
frequencies  reflect  rate  dependence  of  the  hydrodynamic  torsional 
spring  term  on  the  bluff  body. 
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IV.   ANALYSIS 

A.   STATIC  STABILITY 

The  MUS  stability  in  the  pitch  mode  may  be  interpreted  from  Figure 
50;  i.e.,  the  pitching  moment  coefficient  versus  pitch  angle  plot. 
This  plot  was  obtained  only  for  zero  degrees  of  yaw.   At  minus  ninety 
degrees  pitch  angle,  i.e.,  the  ascending  mode,  the  MUS  will  trim 
but  be  statically  unstable.  While  the  pitching  moment  itself  is 
zero  at  o(  =  -90  ,  any  increase  in  pitch  angle  results  in  a  positive 
pitching  moment  which  tends  to  increase  the  pitch  angle  further. 
Similarly,  a  decrease  in  pitch  angle  results  in  a  negative  pitching 
moment.   At  plus  ninety  degrees  pitch  angle  i.e.  the  descending  mode, 
the  MUS  is  again  statically  unstable.   In  the  neighborhood  of  zero 
degrees  pitch  angle,  the  MUS  has  a  positive  pitching  moment  coef- 
ficient and  hence  is  not  trimmed.   Therfore,  in  the  absence  of  the 
cable  constraint,  the  MUS  is  statically  unstable  in  the  pitch  mode. 
The  balsa  wood  dynamic  model  was  suspended  in  a  manner  so  as  to  allow 
a  pitch  degree  of  freedom  for  one  set  of  exploratory  investigations. 
Its  dynamic  behavior  tended  to  confirm  the  preceding  conjectures  since 
the  stable  trim  points  appeared  to  be  weak  and  the  model  was  quite 
prone  to  rotate  in  pitch. 

In  the  vertical  translation  mode  the  MUS  is  statically  unstable. 
Referring  to  Figure  51  f  o{   =0,    y    =  0,  t***0?>  i-n  a  direct  cross  current 
the  MUS  has  a  negative  lift,  which  will  tend  to  cause  the  MUS  to 
descend.   This  in  turn  would  increase  the  effective  pitch  angle, 
and  result  in  a  more  negative  lift. 
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Figure  53  shows  that  the  station  is  statically  stable  in  the  yaw 
mode  (at  zero  pitch)  at  yaw  angles  of  zero  and  180  degrees  and 
statically  unstable  in  yaw  at  90  and  270  degrees  of  yaw. 

Figure  55  shows  that  the  station  is  statically  unstable  in  lateral 
translation  at  yaw  angles  of  zero  and  180  degrees.   If  a  lateral 
translational  velocity  occurs,  an  effective  yaw  angle  is  introduced 
in  the  nature  of  sideslip.  When  this  occurs,  the  side  force  is 
such  that  the  lateral  translation  tends  to  be  further  increased 
(forced) . 

B.   DYNAMIC  STABILITY 

In  the  absence  of  constraints,  the  experiments  indicated  that  the 
MUS  would  be  dynamically  unstable  in  all  translation  modes  and  in 
the  pitch  mode. 

The  yaw  mode  appeared  to  be  dynamically  stable,  as  the  oscillation 
amplitude  either  decreased  in  time  or  remained  constant.   Negative 
damping  ratios  were  not  observed  for  the  range  of  tunnel  speeds  used, 
but  it  is  possible  that  higher  dynamic  pressures  could  have  induced 
negative  damping. 

The  scaling  relations  for  the  yaw  mode  and  the  lateral  translation 
mode  indicate  that  the  natural  frequencies  of  these  modes  are  very 
widely  separated  at  full  scale.   This  would  indicate  that  the  coupling 
tendency  observed  in  the  model  experiments  would  not  occur  for  the 
full  scale  MUS.   Thes  most  important  conclusion  from  scaling 
considerations  is  that  an  oscillation  on  the  model  in  a  wind  tunnel 
on  the  order  of  a  second  per  cycle  period  would  scale  to  the  prototype 
vehicle  in  water  with  a  period  on  the  order  of  20  to  50  minutes  per 
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cycle.   Therefore  it  is  quite  possible  that  unusual  oscillational 
behaviors  in  the  tunnel  would  not  be  perceptible  to  the  crew  enclosed 
in  the  MUS  because  of  the  low  velocities  involved. 

C.   COMMENTS  ON  MUS  CONSTRAINTS 

The  experimental  results  indicated  that  the  proposed  free  ascent 
mode  would  be  practicable.   Considering  an  assumed  stabilized  condition 
where  the  MUS  would  be  floating  freely  in  any  existing  current,  with 
cable  and  pads  jettisoned,  and  in  the  absence  of  an  unequilibrated 
buoyancy  force,  no  hydrodynamic  force  would  be  generated.   If  the 
positive  buoyancy  force  of  approximately  100  lbs.  were  considered, 
the  MUS  would  be  forced  into  vertical  translation  and  hydrodynamic 
forces  would  be  generated.   In  this  orientation,  the  wind  tunnel 
axis  position  of  minus  ninety  degrees  of  pitch  angle  would  be 
applicable  (Fig.  2).   Assuming  that  a  constant  side  current  existed 
throughout  the  water  column,  the  MUS  would  begin  to  rise,  sensing 
only  the  hydrodynamic  forces  corresponding  with  o(  =  -90  . 
However,  in  this  orientation,  both  a  drag  force  and  a  lift  force  could 
be  generated.   If  one  assumes  that  the  MUS  could  achieve  a  vertical 
velocity  of  0.2  knots,  from  the  wind  tunnel  results,  the  resulting 
drag  force  would  be  100  lbs.,  and  the  lift  force  would  be  zero 
(Figs.  49,  51,  52).   Actually,  the  resulting  drag  forces  would  limit 
vertical  ascent  speed  to  0.2  knots,  if  the  100  lb.  buoyancy  force  were 
maintained.   This  speed  limitation  may  be  verified  by  cross  referencing 
Figs.  49  and  52.   The  drag  coefficient  of  0.7  will  result  in  a  drag 
force  of  100  lbs.  at  0.2  knots,  therefore  the  drag  force  will 
equilibrate  the  net  buoyancy  force  at  0.2  knots.   A  check  of  the 
pitching  moment  vs.  pitch  angle  plot  (Fig.  50)  for  this  orientation 
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shows  that  the  MUS  would  be  at  an  unstable  trim  point,  that  is  the 
pitching  moment  would  be  zero  but  the  variation  of  moment  with  flow 

angle  would  not  produce  a  restoring  moment.   As  long  as  the  flow 

I  ° 

vector  remains  constant  at  Q\    =  -90   ,  the  MUS  will  remain  upright. 

If  the  flow  vector  is  disturbed  from  0(  =  -90   or  more  likely,  if 
the  MUS  begins  ascending  from  an  off-vertical  orientation,  the 
resulting  pitching  moment  will  be  in  a  direction  tending  to  over- 
turn the  MUS.   This  pitching  moment  would  be  on  the  order  of  300 

o 
ft-lbs.  maximum  (for  a  0.2  knot  vertical  speed  and  a  20  off -vertical 

orientation).   Considering  the  design  (buoyancy  couple)  stability, 

this  slight  pitching  moment  should  present  no  difficulties  (ref.  1). 

That  is,  the  location  of  the  center  of  buoyancy  above  the  center 

of  gravity  results  in  a  stabilizing  moment. 

Considering  the  case  where  the  MUS  is  tethered  by  the  cable  in 
the  mid-water  column,  if  side  currents  exist  with  a  speed  of  0.3 
knots,  a  pitching  moment  would  be  generated  on  the  order  of  1300 
ft-lbs.  maximum.   This  would  place  a  side  load  on  the  cable  attach- 
ment point,  as  the  cable  would  equilibrate  the  pitching  moment.   In 
addition,  a  yawing  moment  would  exist,  having  the  same  order  of 
magnitude.   The  yawing  moment  would  apply  a  torsional  load  to  the 
cable  or  to  the  cable  attachment  point,  depending  on  the  actual 
design,  however,  the  vehicle  would  tend  to  align  itself  to  the 
stable  trim  point  corresponding  to   1u  =   0  or  180. 

The  oscillatory  behavior  of  the  model  appears  to  be  of  no 
practical  consequence  to  the  full  scale  MUS  and  probably  would  not  be 
noticable  by  the  crew.   Considering  the  proposed  emplacement  procedure, 
since  the  principal  static  equilibrating  force  that  exists  to  counteract 
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drag  forces  (in  side  currents)  is  the  100  lbs.  buoyancy  force,  the 
wind  tunnel  measurements  indicate  that  a  forced  descent  would  result 
with  current  speeds  as  small  as  0.2  knots.   At  this  speed,  the  resulting 
drag  force  would  be  equal  to  the  buoyancy  force,  and  therefore,  static 
equilibration  would  not  be  possible. 


The  restoring  force  F   is  equal  to 

x 


F  =  T  sin  ^ 
x 


Where  T  =  buoyancy  force  =  cable  tension 

and  ^  =  deflection  angle  of  the  cable 

from  the  vertical 

if  T  =  100  lbs. 

for  a  current  speed  of  0.15  knots  and  C   =  0.7 

D  =  60  lbs.  (Fig.  52) 

Sin  ^   =  0.6 

^=37° 

or  in  other  terms,  a  horizontal  translation  of  3,600  ft.  and  a  1,200 

ft.  descent  would  occur  before  static  equilibration  if  6,000  ft.  of 

cable  separated  the  MUS  and  its  anchor.   For  this  reason,  it  is 

recommended  that  a  higher  buoyancy  force  be  considered. 

The  foregoing  discussion  on  forced  descent  is  admittedly  overly 

simplified  and  does  not  take  into  consideration  the  special  effect 

of  catenary  behavior  of  the  cable  itself.   This  topic  would  merit 

individual  study  and  analysis  before  a  complete  discussion  of 

forced  descent  could  be  made.   It  is  also  recognized  that  increasing 

the  buoyancy  force  will  create  a  serious  problem  involving  free 

ascent  terminal  velocity.   In  this  respect,  a  compromise  would  be 

involved  between  stability  considerations* and  ascent  problems. 
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TABLE  I 

VALUES  OF  SOME  PHYSICAL  CONSTANTS 

2 

Characteristic  area  (model) 0.546  ft 

2 

Characteristic  area  (prototype)   .......  1248  ft 

Height  of  model , 1.0  ft 

Height  of  prototype 47.8  ft 

Characteristic  length  (model)   .  .  3.75  in 

Characteristic  length   (prototype)  ...  ,  .  14.93  ft 

Scale  ratio,  A, 0.0209 

vl  -4 

A   4.37  x  10 


h 9.0  x  10 

v5"  -9 

>S 4.0  x  10 

-4   2 
~0,  salt  water 0.17  x  10   ft  /sec 

/    0 

"O,  air 1.58  x  10   ft  /sec 


H 


1.23  x  10"3 

Dimensionless  velocity,  V 443 


Note:   The  characteristic  length  was  arbitrarily  chosen  to  be  the 
average  or  mean  of  the  diameters  of  the  five  principal 
components  which  comprise  the  MUS. 
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TABLE  VI:   RUN  SCHEDULE 


A.   Rigid  Model  Tests,  Model  A  in  Aerolab  Wind  Tunnel 


Dy: 

namic 

est 

Run 

Pr< 

assure, 

Pitch 

Yaw 

Numb 

er 

(psf) 

Angle 

Angle 

Comments 

1 

3.85  to  57.8 
in  increments 

0 

0 

^A 

C     VS 

>  M 

of 

4.0 

2 

ii 

-30 

0 

ii 

3 

it 

+30 

0 

ii 

4 

ii 

0 

30 

ii 

5 

ii 

-30 

30 

it 

6 

ii 

+30 

30 

ii 

7 

ii 

0 

60 

ii 

8 

M 

-30 

60 

it 

9 

II 

+30 

60 

ii 

10 

II 

0 

90 

ii 

11 

II 

-30 

90 

it 

12 

II 

+30 

90 

ii 

13 

II 

0 

120 

ii 

14 

II 

-30 

120 

ii 

15 

II 

+30 

120 

ii 

16 

II 

-30 

150 

ii 

17 

II 

0 

150 

ii 

18 

II 

+30 

150 

ii 

19 

II 

0 

180 

ii 

20 

II 

-30 

180 

ii 

Re 
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Dynamic 

Test 

Run 

Pressure 

• 

Pitch 

Yaw 

Numb 

er 

(psf) 

Angle 

Angle 

Comments 

21 

3.85  to  57.8 
in  increments 

+30 

180 

c>  fi.fi*   «  «& 

of  4.0 

22 

ii 

0 

210 

23 

it 

-30 

210 

24 

ii 

+30 

210 

25 

it 

-30 

240 

26 

ii 

0 

240 

27 

ii 

+30 

240 

28 

ii 

0 

270 

c*  vs  *v 

29 

ii 

-30 

270 

30 

ii 

+30 

270 

31 

ii 

0 

300 

32 

ii 

-30 

300 

33 

ii 

+30 

300 

34 

ii 

0 

330 

35 

ii 

-30 

330 

36 

ii 

+30 

330 

37 

46.1 

-30 

to  +30  in 

0 

co,C  CM  VS  o< 

5 

increments 

. 

>    ■ 

38 

3. 

.85  to  57. 

8 

+90 

0 

c0CLC„  vs  WR. 

in  increments 

of  4.0 

39 


40 


46.1      +60  to  +120  in 
5   increments 

3.85  to  57.8       +60 
in  increments 
Of  4.0 


180 


44 


Test  Run 
Number 

Dynamic 
Pressure, 
(psf) 

Pitch 
Angle 

Yaw 
Angle 

Comments 

41 

3.85  to  57.8 
in  increments 
of  4.0 

+60 

0 

Co 

,CSCM  VS  Nt 

42 

ii 

+60 

30 

ii 

43 

ii 

+60 

210 

H 

44 

34.62  to  46.1 

+90 

30 

ti 

45 

3.85  to  42.3 
in  increments 
of  4.0 

+60 

60 

ii 

46 

ii 

+60 

240 

ii 

47 

30.8  to  38.5 

+90 

60 

it 

48 

3.85  to  42.3 

+60 

90 

ii 

49 

3.85  to  38.5 

+90 

90 

ii 

50 

3.85  to  42.4 
in  increments 
of  4.0 

+60 

270 

^ 

>,Cl,Cm  Ki, 

51 

ii 

+60 

120 

it 

52 

ii 

+60 

300 

ii 

53 

ii 

+60 

150 

H 

54 

ii 

+60 

330 

tt 

55 

46.1 

in 

-120  to  -60 
5   increments 

0 

C>, 

^,CM  vs  « 

56 

3.85  to  46.1 
in  increments 
of  4.0 

-60 

180 

^ 

AA.**! 

57 

ii 

-90 

0 

ii 

58 

ii 

-60 

0 

ii 

59 

ii 

-60 

210 

ii 

«e 


Re 


45 


Dynamic 

Test  Run 

Pressure 

Pitch 

Yaw 

Number 

(psf) 

Angle 

Angle 

60 


61 


62 


63 


64 


65 


66 


67 


68 


69 


70 


71 


46.1 


345  to  015  in 
5   increments 

020  to  045  in 
5   increments 

050  to  075  in 
5   increments 

080  to  105  in 
5   increments 

110  to  135  in 
5   increments 

140  to  165  in 
5   increments 

170  to  195  in 
5   increments 

200  to  225  in 
5   increments 

230  to  255  in 
5   increments 

260  to  285  in 
5   increments 

290  to  315  in 
5   increments 

320  to  340  in 
5   increments 


Comments 
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B.   Dynamic  Model  Tests,  West  Coast  Research  Wind  Tunnel 


Photographs        Time-History 
(Motion  Picture)        Trace 


none  none 

yes  none 

no  yes 

no  yes 

no  yes 

no  yes 

no  yes 

no  yes 

no  yes 

no  yes 

no  yes 

no  yes 

yes 

film  test  yes 

sequency  16  (note  5) 

25A         C        3.0  yes  yes 

film  test 
sequence  17 

26  C        4.0  yes  yes 

film  test 
sequence  18 

27  C        5.0  yes  yes 


Test  Run 

Model 

Dynamic 

No. 

Used 

Pressure, 

Not  Numbered 

0.5  to 

(Note 

1) 

B 

1.0  psf 

M-l  -  M-15 

(Note 

2) 

C 

5.0 

1-9 

(Note 

3) 

C 

5.0 

10,  11 

C 

2.5 

12 

C 

7.5 

13,  14 

C 

2.0 

15,  16 

C 

3.0 

17,  18 

C 

4.0 

19,  20 

C 

5.0 

21,  22 

C 

6.0 

23 

C 

7.0 

24,  25 

C 

8.0 

24A 

(Note 

4) 

C 

2.0 
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Test  Run 

Model 

Dynamic 

Photographs 

Time-History 

No. 

Used 

Pressure, 

(Motion  Picture) 

Trace 

28 

C 

5.0 

yes 
film  test 
sequence  19 

yes 

29 

C 

6.0 

yes 
film  test 
sequence  20 

yes 

30 

C 

7.0 

yes 

film  test 

sequence  21 

yes 

31 

C 

7.0 

yes 
film  test 
sequence  22 

yes 

32 

C 

8.0 

yes 
f  i.lm:  test 
l sequence  23 

yes 

33  C        8.0  yes  yes 

no  yes 

no  yes 

no  yes 

no  yes 

no  yes 

no  yes 

no  yes 

no  yes 

no  yes 

no  yes 

no  yes 


L-l,  L-2, 

L-3 

D 

2.0 

L-4,  L-5 

D 

3.0 

L-6,  L-7 

D 

4.0 

L-8,  L-9 

D 

5.0 

L-10,  L-ll 

L-12 

D 

6.0 

L-13,  L-14 

D 

7.0 

L-15,  L-16 

L-17,  L-18 

D 

8,0 

LL-1,  LL-2 

E 

2.0 

LL-3,  LL-4 

E 

3.0 

LL-5,  LL-6 

E 

4.0 

LL-7,  LL-8 

LL-9 

E 

5.0 
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Test  Run     Model  Dynamic  Photographs  Time-History 

No.       Used  Pressure,  (Motion  Picture)         Trace 

LL-10,  LL-11   E        6.0  no               yes 

LL-12,  LL-13 

LL-14,  LL-15   E        7.0  no               yes 

LL-16,  LL-17   E        8.0  no               yes 


Notes: 

1.  Various  runs  conducted  in  both  the  Aerolab  tunnel  and  the 
West  Coast  Research  tunnel.  Model  mounted  on  wire  in  order 
to  observe  both  the  yaw  mode  and  the  pitch  mode. 

2.  Local  tie  down  constraints  not  used;   these  runs  correspond 
to  Part  One  of  the  movie  film.   Film  test  sequences  1 
through  15;  film  speed,  16  frames  per  second. 

3.  Local  tie  down  constraints  were  used  in  this  and  subsequent 
runs  listed. 

4.  Runs  24A  through  33  correspond  to  Part  Two  of  the  movie 
film.   Film  test  sequences  16  through  23;  film  speed, 
16  frames  per  second. 

5.  Test  run  numbers  (24A  through  33,  L-l  through  L-18,  LL-1 
through  LL-17)  correspond  with  notation  used  on  the 
time-history  traces. 
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1.  MODEL      SHOWN      AT     ZERO      PITCH       AND      ZERO    YAW, 

2.  FORCES      AND     MOMENTS       ARE      RELATIVE       TO 
WIND       AXIS      SYSTEM     (AS      ABOVE). 


FIG.  2 


COORDINATE      AXIS      SYSTEM 
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FIG.    3  YAW     ANGLE      DEFINITION 
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FIG. 4  PITCH     ANGLE      DEFINITION 
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FIG 

5 

EXTERIOR 

DIMENSIONS 

OF 

THE 

MUS 

MODEL 

54 


3] 
Q 
CD 


C 


c 

o 

o 

oc 

m  en 

r-m 

D 


I 

m 


55 


FIG.    7 

WEIGHT 

DISTRIBUTION 

IN     THE 

DYNAMIC     MODEL 

WEIGHTS      A.B.C  ,D,E  ,F, 
AND     1      -     1.0   in.  dia 
G  ,H   -    2.0  in.    dia 


0.43        f 

"  f    -   I 
I         0.9  G 


WEIGHT     B 


Jf    _    f_BASE 
LINE 
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COMPONENT 

WEIGHT 

DIA. 

HEIGHT 

WEIGHT 

A 

0.30 

lbs 

1.0 

Q96 

WEIGHT 

B 

0.67 

lbs 

1.0 

207 

WEIGHT 

C 

0.67 

lbs 

1.0 

2.07 

WEIGHT 

D 

0.67 

lbs 

1.0 

2.07 

WEIGHT 

E 

0.67 

lbs 

1.0 

2.07 

WEIGHT 

F 

0.67 

lbs 

1.0 

2.04 

WEIGHT 

G 

1.45 

lbs 

2.0 

1.12 

WEIGHT 

H 

1.00 

lbs 

2.0 

0.77 

WEIGHT 

1 

0.31 

lbs 

1.0 

0.96 

HOLDING 
TANK 

0.20 

lbs 

REACTOR 
CYLINTER 

0.74 

lbs 

HABITAT 
CYXINDEfi 

O.Pr 

lbs 

MAIN    BAt. 
TANK 

.AST 

1.25 

lbs 

OBSERVATION 
SPHERE 

0.05 

lbs 

FIG.    8 

WEIGHT       DISTRIBUTION 
IN    THE     DYNAMIC     MODEL 
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FIG 

9 

FASTENING 

METHOD 

FOR 

THE 

DYNAMIC 

MUS 

MODEL 
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Fig.  11  -  Model  orientation:  yaw  zero  degrees,  pitch 
zero  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


. 


Fig.  12  -  Model  orientation:   yaw  zero  degrees,  pitch 
zero  degrees;  viewed  from  upstream  side 
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Fig.  13  -  Model  orientation:   yaw  zero  degrees,  pitch 
+30  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


Fig.  14  -  Model  orientation:   yaw  zero  degrees,  pitch 
+30  degrees;  viewed  from  upstream  side 
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Fig.  15  -  Model  orientation:   yaw  zero  degrees,  pitch 
-30  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


Fig.  16  -  Model  orientation:   yaw  zero  degrees,  pitch 
-30  degrees;  viewed  from  upstream  side 
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Fig.  17  -  Model  orientation:   yaw  30  degrees,  pitch 
zero  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


Fig.  18  -  Model  orientation:   yaw  30  degrees,  pitch 
zero  degrees;  viewed  from  upstream  side 
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Fig.  19  -  Model  orientation:   yaw  30  degrees,  pitch 
+30  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


Fig.  20  -  Model  orientation:   yaw  30  degrees,  pitch 
+30  degrees;  viewed  from  upstream  side 
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Fig.  21  -  Model  orientation:   yaw  30  degrees;  pitch 
-30  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


Fig.  22  -  Model  orientation:   yaw  30  degrees,  pitch 
-30  degrees;  viewed  from  upstream  side 
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Fig.  23  -  Model  orientation:   yaw  60  degrees,  pitch 
zero  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


Fig.  24  -  Model  orientation:   yaw  60  degrees,  pitch 
zero  degrees;  viewed  from  upstream  side 
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Fig.  25  -  Model  orientation:   yaw  60  degrees,  pitch 
+30  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


Fig.  26  -  Model  orientation:   yaw  60  degrees,  pitch 
+30  degrees;  viewed  from  upstream  side 
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Fig.  27  -  Model  orientation:  yaw  60  degrees;  pitch 
-30  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


Fig.  28  -  Model  orientation:  yaw  60  degrees,  pitch 
-30  degrees;  viewed  from  upstream  side 
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Fig.  29  -  Model  orientation:  yaw  90  degrees,  pitch 
zero  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


Fig.  30  -  Model  orientation;  yaw  90  degrees,  pitch 
zero  degrees;  viewed  from  upstream  side 
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Fig.  31  -  Model  orientation:   yaw  90  degrees,  pitch 
+30  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


Fig.  32  -  Model  orientation:   yaw  90  degrees,  pitch 
+30  degrees;  viewed  from  upstream  side 
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Fig0  33  -  Model  orientation:   yaw  90  degrees,  pitch 
-30  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


Fig.  34  -  Model  orientation:   yaw  90  degrees,  pitch 
-30  degrees;  viewed  from  upstream  side 
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Fig.  35  -  Model  orientation:   yaw  zero  degrees,  pitch 
+60  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


I 


Fig.  36  -  Model  orientation:   yaw  zero  degrees,  pitch 
+60  degrees;  viewed  from  upstream  side 
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Fig.  37  -  Model  orientation:   yaw  zero  degrees,  pitch 
+90  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


Fig.  38  -  Model  orientation:   yaw  zero  degrees,  pitch 
+90  degrees;  viewed  from  upstream  side 
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Fig.  39  -  Model  orientation:   yaw  zero  degrees,  pitch 
-60  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


Fig.  40  -  Model  orientation:   yaw  zero  degrees,  pitch 
-60  degrees;  viewed  from  upstream  side 
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Fig.  41  -  Model  orientation:   yaw  zero  degrees,  pitch 
-90  degrees;  viewed  from  left  side  of  model;  airflow 
left  to  right 


Fig.  42  -  Model  orientation:   yaw  zero  degrees,  pitch 
-90  degrees;  viewed  from  upstream  side 
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Fig.  43  -  Support  arrangement  used  for  yaw  variation 
in  five  degree  increments;  model  at  zero  degrees 
pitch;  viewed  from  side;  flow  left  to  right 


Fig.  44  -  Same  as  43,  but  view  from  upstream  side 
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Fig.  56  -  Photograph  of  the  West  Coast  Research  48  in. 
x  60  in.  wind  tunnel,  showing  support  beams  for  wire 
attachment  exterior  to  the  test  section 


Fig.  57  -  The  dynamic  model  in  the  test  section, 
side  view 
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Fig.  58  -  Photograph  of  the  dynamic  model  in  the 
test  section,  view  from  downstream 


Fig.  59  -  Close-up  of  the  dynamic  model 


Fig.  60  -  Photograph  of  dynamic  model  instrumentation 
hardware  and  bearings 
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Fig.  61  -  Photograph  of  the  modified  potentiometer 
used  in  dynamic  model  instrumentation 
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FIG.  67        TIME    HISTORY    OF    FREE 
VIBRATION    IN     YAW 
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